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TECHNICAL NOTE 2206 


GRAPHICAL METHOD FOR OBTAINING FLOW FIELD IN TWO-DIMENSIONAL 
SUPERSONIC STREAM TO WHICH HEAT 13 ADDED 
By I. Irving Pinkel and John S. Serafini 


SUMMARY 

A graphical method for obtaining the two-dimensional supersonic 
flow with moderate heat addition is presented. The basic equations 
were derived for flows in which the total temperature has a continuous 
variation with displacement in the direction of the flow. For flows 
at constant total temperatures before heat addition and with total 
temperature varying only in the direction of the streamlines, simpli- 
fied forms of the basic equations were obtained from which working 
charts were plotted to facilitate construction of the flow. The 
working charts cover the range of Mach number from 1.3 to 5.0. The 
application of the method and the charts are illustrated by a simple 
example, giving the flow around a curved surface with addition of 
heat in an adjacent finite layer where the initial total temperature 
of the flow is the same on all streamlines and lines of constant 
total temperature in the heated portion of the flow lie approximately 
normal to the streamlines. 


INTRODUCTION 

Current interest in supersonic flow includes flow through gas 
turbines and in jet-engine exhausts in which heat release by contin- 
uing combustion of fuel introduced in an upstream combustion chamber 
modifies the flow. The general problem of the effect of heat addition 
on supersonic flows is of growing interest in the supersonic-flight 
and missile fields. The addition of heat to a gas flowing in a 
constant-area duct with a steady supersonic velocity has been investi- 
gated (reference l) . A theory for one -dimensional steady compressible 
fluid flow in ducts with heat addition is presented in reference 2 . 

An extension of the method of characteristics (for example, see refer- 
ence 3) for the graphical construction of irrotational supersonic 
flows was developed at the NACA Lewis laboratory and is presented 
herein for obtaining the two-dimensional supersonic flow to which 
heat is added. The method is restricted to flows in which the heat 
transfer between adjacent stream tubes can be neglected and entropy 
changes in the flow are due only to the addition of heat (shock- 
free flow) . 
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The method of characteristics provides a construction of the 
irrotational flow field subdivided into zones in each of which the 
Mach number, and therefore all properties of the flow, are constant. 
When heat is added to the flow, further division of the flow is 
required to account for the variation of total temperature and total 
pressure throughout the field. In the rotational flow that results 
from the addition of heat, compression and expansion waves arise in 
the flow to maintain continuity of static pressure across adjacent 
streamlines. Because of the changes in total pressure and tempera- 
ture and the generation of waves within, the flow streamlines and 
heat-addition reference lines must be plotted with the waves. The 
complexity of the graphical construction of the flow that results 
from heat addition and wave generation can be reduced by methods 
presented herein. 

The equations developed as the basis of the proposed method ap- 
ply to any mode of heat addition that gives a continuous variation of 
total temperature with displacement in the direction of the flow (but 
not necessarily continuous variation normal to the streamlines) . For 
flows at constant total temperature before heat addition and with 
constant total temperature along lines approximately normal to the 
streamlines, simplified forms of the basic equations are obtained 
from which working charts have been plotted to facilitate the con- 
struction of the flow. 

As an illustrative example, application of the proposed method 
is made to the problem of obtaining the heated supersonic flow field 
adjacent to a curved wall. The initial total temperature of the flow 
is the same on all streamlines and lines of constant total tempera- 
ture in the heated portion of the flow lie approximately normal to 
the streamlines. 


SYMBOIS 


The following symbols are used in this report: 


c. 


P 


specific heat at constant pressure, herein assumed to be 
constant 


c. 


v 


specific heat at constant volume 


M 


Mach number 


N square of Mach number 
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P total pressure 

p static pressure 

Q heat content of flow based on total temperature, Q = CpT 

S gas constant 

T total temperature 

t static temperature 

v velocity 

a angle through which supersonic stream is turned to accelerate 

from M = 1 to M > 1 

Aa wave strength expressed in terms of deflection produced in 
supersonic stream (positive value of Aa is that which 
increases static pressure) 

p Mach angle, siri"^ l/M 

7 ratio of specific heats, Cp/cy 

9 angle between local velocity and direction of initial flow 

p density 

Subscripts: 

0 refers to initial conditions before addition of heat 

1,2,3,... designate particular sector of flow field bounded by 
waves, by heat-addition lines, or by both 


METHOD 

Basic Considerations 

For the purposes of the method presented, the addition of heat 
to the flow in a stream tube is considered a two-step process. The 
stream tube is divided into short intervals along its length, over 
each of which a small change in heat content of the flow occurs. All 
the heat addition in each interval is considered to take place along 
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a line at the middle of the interval. The heat addition is assumed 
to take place at constant stream-tube area. Following the heat addi- 
tion, the streamlines change direction, with a corresponding change 
in stream-tube area to adjust the local static pressure as required 
by the adjacent flow. 


Governing Equations 

The change in static pressure with a heat addition AQ at con- 
stant stream-tube area is obtained by application of the conservation- 
of -momentum and the continuity equations. For the constant-area case, 
these equations have the form 

p + pv 2 = K x (1) 

pv = K 2 (2) 

where Kq and K 2 are constants. From the equation of state for 
perfect gases, 


p = pBt 


and the equation for Mach number, 

M = -L_ 

(3) 

(4) 


equations (l) and (2) become, respectively, 

p(l+yN) = Kq 

P+. ^ 


where N = M 2 . Differentiation of equations (3) and (4), written in 
the approximate-difference form, and substitution of 

T = t(l + N | yield, on the elimination of AN/N and Ap/p, 

successively, 


P " Tn^T) t 



( 5 ) 
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where 


AN _ ( 1+yN )(l + 2 K ) AT _ D AQ 

N (N-l) T “ N Q 


7N 


C = 


( x + ¥ N ) 


Tn^TT 


N(1+ 7 N) 


D = - 


» {* ¥ H ) 

N-l 


(6) 


(7a) 


(7b) 


Q = c p T (7c) 

^ = AZ (74) 

Q T 


Adjustment of the pressure following heat addition, as required by 
the adjacent flow, is accomplished by a wave of strength Aa gen- 
erated in the flow at the point of adjustment shown in figure 1. 

If Ap is the required pressure change across the wave, the value 
is given by the simplified relation suitable for waves of small 
strength, 


Ap = Aa = pA Aa 

-v/N-1 


( 8 ) 


where 


A = 4^ (8a) 

yn 

and a positive value of Aa is that which produces an increase in 
static pressure. The sura of static-pressure changes resulting from 
heat addition at constant area, followed by pressure adjustment, is 
obtained from equations (5) and (8): 


AP ’4+ 2 T I, ) m 

p " (N-l) Q 


+ ■ ? — 4a = C ^ + A Aa 
v/N-1 ft 


( 9 ) 
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The Mach number change in the flow crossing the wave may be 
obtained, by differentiating the well-known Prandtl-Meyer expression 



to obtain the expression for a small finite change in a 


or 


-Aa 


an yin 

2N 1 + 111 N 


AN 

N 


2 ( lt 2 r”) 

VN 7 ! 


Aa 


(10a) 


(10b) 


The total change in AN for the heat addition at constant area and 
subsequent pressure adjustment is obtained from equations (6) and 
(10b) as 


AN 
N ~ 


(UrH)(l + 2|i ») aQ 
N-l Q 


2 ( 1 + 2 t m ) 

Vn 7 ! 


(10c) 


Application of Governing Equations 

Equations are now developed for four distinct cases of flow 
adjustment that arise in the solution of heated supersonic flows by 
the method of this report. To each case of flow adjustment corre- 
sponds a generated wave pattern that can be computed by application 
of equations (5) to (10b). Because the method of solution applies 
to any mode of heat addition that provides a continuous variation of 
temperature along a streamline , the heat content per unit mass of 
fluid Q may vary discontinuously with displacement across stream- 
lines. For the special case in which the total temperature is con- 
stant on all streamlines upstream of the heat-addition region and 
lines of constant total temperature lie normal to the streamlines in 
the heat-addition zone, simplified expressions suitable for computing 
working charts are obtained for the waves generated in the four cases 
of flow adjustment. 
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The marked effect of small quantities of heat added to the 
flow , representing a change of several percent in total tempera- 
ture or heat content Q, on the local properties of the flow 
requires that heat-addition reference lines be drawn in the flow 
to devide the flow into zones in which a change in heat content Q 
of only 2 or 3 percent occurs. These reference lines, drawn normal 
to the local streamlines for convenience, represent the interval 
center lines discussed in the previous section, at which all the 
heat added in the interval is concentrated. 

Case I: wave generation at boundary between heated and unheated 

flow . - The first case considered is schematically illustrated in 
figure 1 , which shows the flow adjustment that occurs along a stream- 
line separating heated and unheated portions of the flow. The pres- 
sure equilibrium that exists across the bounding streamline before 
heat addition is restored after heat addition by a deflection of the 
flow accompanied by the generation of two waves of equal magnitude 
but opposite sign. Subscripts used in the following derivations 
refer to flow regions numbered in the figure. From equation (5), 
the pressure P 3 in the heated stream tube, after a head addition 
of AQ, is 


where (AQ )2 and Qg are associated with the flow in region 2. If 
Aap is the turn in the bounding streamline required to restore pres- 
sure equilibrium in the regions 4 and 5, then from equation ( 8 ), 


If Aaq is the wave strength associated with the flow having the 
value of N = Nq, then the wave separating regions 3 and 4 has the 
strength -Aaq. With the condition that the static pressures on 



(Ha) 



(Hb) 


3 



(11c) 
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both, sides of the bounding streamline are equal after a heat addition 
of (AQ) 2 , p 4 = p^, and with the substitution of equation (11a) in 

equation (lib), the following expression is obtained: 


Pl + 51 (Jk\ = P 2 + ^ 


rN 


i 1 : ¥ M ) 


N-l 


(AQ); 

Q? 


P2 


f 




l 1 _L 

7 n^ 

") 

(W) 2 _ 



N-l 

^2 




2 





Ao^ 


(lid) 


With the substitution in equation (lid) of N3 in terms of N2 
by means of equation (6), and from the relation p-^ = p 2 , the solution 
of the equation (lid) for Aa^ in terms of known quantities is 


C 2 


Aaq = 


(AQ); 

Q 2 


A ! + 


(AQ), 


1 + C c 


7F. 


( 12 ) 


Vvi 


where A, C, and D are defined by equations (8a), (7a), and (7b), 
respectively, and 


F = N + D 


AQ 

Q 


(12a) 


The expression for Aap (equation (12)) can be written as a function 
of Np/N2 and N2 that is convenient for presentation in chart form. 
The chart and its uses are presented when the complete graphical 
method is considered. 

The pressure change p 4 -p 2 l 8 obtained by means of equations (11c) 
and (12) as 
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P4-P2 = P5-P1 = Pi < 


r 

) 

A1C2 

<a«) 2 
- Q 2 - 


Ai + 

, „ (AQ)a 

1 + Co 7T 

£ Q2 

7*2 

VF2-1 


y 


(13) 


Case II; Interaction between wave and heat -add It Ion zone. - In 
the portion of the flow schematically depicted in figure 2. a wave 
crosses a heat -addition zone at point P. The flow above the stream- 
line through P crosses the incident wave before encountering the 
heat addition zone. The flow in region 2 below the streamline through 
P crosses the heat-addition line first. In order to equalize the 
pressure in regions 4 and 5, the incident wave, on transmission 
through the heat-addition zone at P, results in a modification of 
the incident wave to a strength Aa^ and a generation of a wave of 
strength Aa^. The pressure in region 7 is given (from equation (8)) 
as 


P7 = Pi + Pi 



Aa 


(14a) 


The corresponding change in Mach number may be obtained by use of 
equation (10b). The pressure in region 6 is obtained by means of 
equation (5) in terms of the flow quantities in region 7; 


7N 1 


p = p + p 7 


+ ¥ ■) 


N-l 


(AQ) X 

“oT 


(14b) 


where AQq and Qi are associated with the flow from region 1. 

In order for the flows in regions 4 and 5 to have the same 
direction, the strength of the wave Aaj, arising at point P, is 
determined as 


Aaj = Aa - Aajj- 


(15) 


The pressure in region 5 is given in terms of the flow quan- 
tities in region 6 as 
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P5 = P6 



(Aa-Aa^.) 


(14c) 


By substituting successively for the pressure p 6 and P7 ; the 
pressure in region 5 may be obtained in the following form: 



(I4d) 

When equations (10b) and (6) are used to obtain the values of N at 
regions 7 and 6, the equation (l4d) takes the form 


P5 = Pi ^1 + A^Aa^ (l +' 


7 ( N+BAa) |l + III ( N+BAalh ( AQ ) ± 




N+BAa-1 


Ql 



(AQ)- 


7 N+BAa - 


( t [ 

(N+BAa) (1+yN+yBAa) [l4-2^1(N+BAajJ q - 


N+BAa -1 


N+BAa-1 - 


(N+BAa) (l+yN+yBAa) |l+Ill( N+BAa) 


(Aa-Aofc) 


N+BAa - 1 


(I4e) 


where 


2N 


B = - 


( 1 < 2 r') 

yin 


(16) 


Equation (16) may be more simply written, by means of the substitution 


E = N + B Aa 


(I4f) 
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P 5 = p 2 (1 + AjAa.) 


n + 

7E (l + e) 

(AQ)i 1 


E-l 

1 Ql J 


1 + 



(I4g) 


In a similar manner, the pressure in region 4 may he obtained, as 


P 4 = P2 


h. 


(AQ)2 

~^T 


+ 7 


1 + C, 


(AQ); 


Fc 


q 2JVf^I 


Aa k > 


(I4h) 


where C and. F are functions of N as defined, in equations (7a) 
and (12a), respectively. With the use of pq = P 2 and P 5 = p 4 , 
equations (I4g) and (I4h) may be solved for Aa k as a function of 
the known quantities Ny, N 2 , (AQ)i/Qi, (AQ)2/Q2> and Aa: 


Z n -1 


V. 


(17a) 
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where 




f (18) 


Where the initial total temperature is constant and the heat addi- 
tion is considered constant along lines normal to the local stream- 
lines , = Qg = Q and (AQ)q = (AQ) 2 = AQ. Therefore, with 

Ni = N 2 = N, equation (17a) is written in terms of N, AQ/Q, and 
Aa as 


A at 


Z-l + YZ Aa - C ^ 

Q_ 

YZ + W 


(17b) 


where Y, Z, and W are evaluated as given by equations (18), with 
Ql = Q 2 = Q, (AQ)i = (AQ )2 = AQ, and Nq = N 2 = N. The pressure 
change may be written as 


P5-P1 

Pi 


P4-P1 „ (^2 .. 

= C2 — x — + w 2 


Pi 


«2 


Zq-1 + YqZq Aa - C 2 


(AQ)2 

Q2 


YqZq + W 2 


( 19 ) 
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Case III: wave propagation across adjacent stream tubes of 

differing total pressure. - A wave passing through adjacent stream 
tubes of differing total pressure produces unequal changes of static 
pressure in the stream tubes. In order to maintain continuity of 
static pressure and flow direction downstream of the wave, the wave 
is modified in strength and another wave is generated at the stream- 
lines separating two adjacent stream tubes of differing total pres- 
sure. The situation schematically illustrated in figure 3 is con- 
sidered: A wave of strength Aa is shown crossing the boundary 

between regions 1 and 2 in which the flow is parallel and at the same 
static pressure upstream of the wave, but differs In Mach number. If 
Aa-t is the strength of the wave transmitted through the bounding 

streamline and Aa r the strength of the wave generated at the bounding 
streamline, then 


The generated wave Acij. will produce a change of pressure that is 
given by 


P 5 = Pi + 



(20a) 



2 


(20b) 



(20c) 


3 


In order for the flow in adjacent stream tubes to have the same 
direction. 


AOj, = Aa - AOj. 


( 21 ) 


The requirement that the pressures in regions 4 and 5 be equal pro' 
vides, by means of equations (10b), (20), and (21), 
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(I 4 A 0 Aa,) + (l+Ao Act ■) | — ===\ (Acl-Aclo.) ( 22 ) 

W&i U 


Pl(l+Ai Aoj.) = P2 
Inasmuch as = pg, solving for Aat/Aa gives 


Aat 

Aa 


Aa * QiA <UAz aa) 

A 1 + f-2=z\ ( 1+A 2 4a ) 


(23a) 


Where the strength of the incident wave Aa can he assumed to he 
small (Ida | sl.00°), calculations have shown that the pressure and 
the Mach number changes from region 2 to region 4 can he considered 
to result from a combined wave of strength Aa + Aa^. acting on the 
flow in region 2. Equation (23a) then is of the simpler form 


Aat 

Aa 


Al 

a 2 


» 1 - 


— i + 1 


ACty 

Aa 


(23b) 


Case IV: heat addition across adjacent stream tubes of dif- 

fering total pressure. - The difference in the change In static pres- 
sure produced by additions of heat to adjacent stream tubes of 
differing total pressure is compensated by waves of equal intensity 
and opposite sign generated at the streamline separating the two 
stream tubes, as illustrated in figure 4. If (AQ)l/Ql And 
(AQig/Qg are values associated with the flow from the regions 1 

and 2, respectively, the pressure in the regions 6 and 3 are obtained 
from equation (5) as 



(24a) 


(24b) 


P3 = *2 + P 2 


2 
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If Aaq is the wave associated with the flow having the value of 
N = Ni, then the wave separating regions 3 and 4 has the strength 
-Aai. The pressures in the regions 5 and 4 in terms of the flow 
quantities in the regions 6 and 3, respectively, are 


P 5 = P 6 + P 6 



Aaq 


(24c) 


P 4 = P 3 



Aaj 


(24d) 


Because p 3 = P 4 , the substitution of equations (24a) and (24b) in 
equations (24c) and (24d), respectively, yields 


Pi 






(AQ)2 7*2 

Q 2_J 



(25) 


where C and F are functions of N previously defined (equa- 
tions (7a) and (12a), respectively). Because pi = p2, equation (25;, 
solved for Aaq, yields 


Aaq 



(26a) 


If the initial total temperature is constant and the heat addition 
is constant along lines normal to the local streamlines, 

(AQ)i = (AQ )2 = AQ and Q-^ = Qg = Q. Then equation (26a) becomes 
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Aai 


(Ci 


c 2 ) 


AQ 

Q 


1 + C 5 


AQ\ ? f 2 

Q 's/Fpi 


+ 


a 


1 + c 


AQ 

1 T 


) 


VFpi 


(26b) 


Working charts. - By means of equations (5), (6), (12), (17b), 

(23b), and (26b), and with y = 1.4, charts were obtained to facili- 
tate the graphical construction of heated supersonic flows where the 
initial total temperature is the same on all streamlines and the 
lines of constant total temperature in the heated portion of the flow 
are normal to the local streamlines. These conditions are approxi- 
mated in ducts and channels having an initially uniform supersonic 
flow in which heat release by chemical reaction is occurring in the 
body of the flow. 

The exact differential form of equation (6), giving the change 
in N with a heat addition of AQ/Q at constant area, has been 
integrated and the results are given in figure 5 for AQ/Q of 0.02 
and 0.03 with AN = Nq-N2 as the ordinate and Nq as the abscissa. - j 

The static- pressure change and the total-pressure ratio across a heat 
addition have been calculated by integrating the exact differential 
form of equation (5) and using the values of N 2 obtained from 
figure 5 for a given AQ/Q. The variation of the static-pressure 
change and the total -pressure ratio across a heat-addition zone with 
Nq, the value of N preceding heat addition are given in figures 6 
and 7, respectively. 


The value of Aaq of case I as a function of N 2 , with Nq/N 2 
as a parameter, is presented in figures 8(a) and 8(b) by means of 
equation (12) for AQ/Q of 0.02 and 0.03, respectively. The value 
of Aaj in case II as a function of Nq with various values of Aa 
is presented in figure 9 by means of equation (17b) for AQ/Q of 
0.02 and 0.03 for both expansion and compression waves. The ratios 
Aoq/Aa and Ao^/Aa of case III are calculated from equation (23b) 
and presented in figure 10 as functions of Aq/Ag. (Equation (23b), 
which is an approximate form of equation (23a), gives sufficient 
accuracy in the range N = 1.5 to N = 25.0 for small values of Aa. ) 

For | Aa| ~ 1.00°, the use of the chart of figure 10 gives values of 
the modified and generated waves, which have less than 0.5-percent 
deviation from the values Obtained from the mere exact equation (23a). 

The variation of A (equation (8a)) with N is shown in | 

figure 11. The value of Aaq of case IV, calculated from equa- 
tion (26b), is presented in figures 12(a) to (I2d) as a function of 
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Ni with various values of W 2 . Figures 12(a) and 12(h) give the 
values for Aa^ of case IV with AQ/Q = 0.02, and figures 12(c) 
and 12(d) give the values for Aa^ of case IV with AQ/Q = 0.03. 
For a given AQ/Q, the curve of Aa,]_ against Nj_ is coincident 
for N 2 = 3.0 and N 2 = 4.0 and very similar for 3.0 < Ng < 4.0. 
Figures 13 and 14 show plots of N against a, the Prandtl-Meyer 
turning angle, and N against p, the Mach angle, respectively, 
with expanded scales that are convenient for graphical plotting. 


COMPLETE GRAPHICAL METHOD 

The initial layout of the graphical construction for heated 
supersonic flows begins as for an irrotational supersonic flow. 

Curved walls hounding the flow are approximated by chord lines of 
small angles of inclination relative to adjacent chord lines. If 
the initial upstream flow is rotational, the flow is divided into 
stream tubes in each of which the Mach number is considered to have 
the value on the streamline through the middle of the stream tube . 

The plot of the flow proceeds by successive application of the four 
cases for which equations were developed in the preceding section of 
the report. It is assumed that in any problem the variation in total 
temperature through the heated portion of the flow is sufficiently well 
known for the accuracy desired in the construction of the flow. As the 
flow construction proceeds, heat -addition reference lines are drawn 
normal to the local streamlines. For flows in which the lines of con- 
stant total temperature also lie normal to the streamlines, the heat- 
addition reference lines are also lines of constant total temperature. 
In the interval between consecutive reference lines, the fractional 
heat addition to the flow AQ/Q should be equal to or less than 0.03 
for values of N from 1.5 to 16.0 in order for the equations repre- 
senting the effect of heat addition on the flow to apply with satis- 
factory accuracy. For values of N greater than 16.0, smaller values 
of AQ/Q are to be used, (for example, for N of 25.0, a satisfactory 
value for AQ/Q is 0.02). When the total temperature is not constant 
along reference lines, the flow is divided into stream tubes, in each 
of which the total temperature is assumed constant and equal to the 
average of the values at the bounding streamlines. A sufficient number 
of stream tubes are chosen to approximate the variation of total tem- 
perature normal to the streamlines . 

The application of the flow-construction method considered is 
essentially the same for flows in which the total temperature is 
constant or varies along lines normal to the streamlines. In order 
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to illustrate the application of the working charts, an example is 
given in which the total temperature is constant on lines normal to 
the streamlines. 

This example, shown in figure 15, presents the contruction of 
a heated supersonic flow adjacent to a curved wall. The flow, 
initially uniform at Mo = 4.0, is hounded by a single wall curved 
through 4° on a circular arc between points A and B and the 
unheated flow extending to infinity. The streamline bounding the 
heated portion of the flow goes through points E and F. A quan- 
tity of heat AQ/Q = 0.03 is added at regularly spaced intervals 
corresponding to 1° turn of the bounding wall. In any interval n, 

Q = Qo(l + AQ/Q) n , where Q q is the value of Q before heat addi- 
tion. Waves arising or reflected from curved walls bounding the flow 
are handled in the manner employed in two dimensional irrotational 
flows. 

The curved wall is approximated by chord lines of increasing 
inclination to the initial flow direction in increments of 1°. The 
first wave is generated at the wall at point C with the strength 
Aa = -1.0°. The corresponding Mach number change and Mach angle 3 
imposed on the flow by the wave are obtained from figures 13 and 14, 
respectively, which give, for Nq = 16.0 and Aa = -1.0°, N2 = 16.62 
and 3 = 14.48°. The value of N and the inclination of the flow to 
the initial flow direction are recorded in figure 15 as N/0, which 
in this case is 16.62/-1.000 0 . Negative values of 0 indicate a 
clockwise turn of the flow where the initial direction of the flow is 
from left to right. At the intersection of the wave and the first 
heat-addition reference line at D (fig. 15), the situation classified 
as case II occurs. The generated and transmitted waves sire obtained 
from figure 9(c) and equation (15). With Aa = -1.00°, Nq = 16.00, 
and AQ/Q = 0.03, Aaj is equal to 0.034° and Aa^ = -1.034°. The 
change in the value of N across the heat-addition reference line 
is obtained from figure 5, with subscripts 1 and 2 representing the 
upstream and downstream flow, respectively. The change in the value 
of N across each wave is obtained from figure 13. The values for 
N and 9 are shown in the appropriate sections of the flow (fig. 15) 
and Aaj and Aat are recorded on the corresponding waves. The 
streamline through point D (fig. 15) is turned through an angle Aa^. 

The flows above and below this streamline have different total pres- 
sures and the streamline DH is drawn in the plot of the flow to 
demark this difference. 
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Along the streamline bounding the heated and unheated flow, two 
waves of equal magnitude but opposite sign are generated at the 
point of intersection E (fig. 15), with a heat-addition reference 
line corresponding to the situation classified as case I. The value 
of Aaq is obtained from figure 8(b) for Nq/N 2 = 16/16 = 1.00, as 
Aaq = 0.899. The change in N of the unheated flow is obtained 
from figure 13 for N = 16 and Aa = Aai = 0.899. The change in N 
of the heated flow across the heat-addition zone is obtained from 
figure 5 for Nq = 16 and AQ/Q = 0.03, as N = 13.33, and across 
the wave Aa = -Aa^ = -0.899 from figure 13 for = 13.33 and 
Aa = -0.899 as N = 13.78. The flow zones involved are labeled with 
the appropriate values of N and 9 and the strength of the waves 
are indicated as well. It should be observed that a positive (com- 
pression) wave extending from the- upper left to the lower right will 
produce a clockwise turning of the flow, as does a negative (expan- 
sion) wave extending from the lower left to the upper right (for 
example, point D). Likewise, a counterclockwise turn of the flow is 
produced by both a positive wave extending from the lower left to 
the upper right and a negative wave from the upper left to the lower 
right (for example, point E). 

At point F in figure 15, a wave is incident on a streamline 
separating regions of differing total pressure, which corresponds to 
case III. An enlarged figure of the zone in the neighborhood of 
point F is shown in figure 16(a). The required values of Aq and 
A 2 are obtained from figure 11 for Nq = 14.55 and N 2 = 10.59, 
giving Aq = 5.534, A 2 = 4.788, and Aq/A 2 = 1.156. By using fig- 
ure 10, the corresponding values of the transmitted and reflected 
waves Aat and Aa r (fig. 3) may be obtained as -0.988° and 
-0.077°, respectively, for Aa = -1.065°. The change in N across 
the various waves is obtained from figure 13 for the corresponding 
initial values of N. The Mach angle between the streamline and 
the waves is obtained from figure 14 as the angle p. It should be 
noted that whenever Aq and A 2 are very nearly equal, the value 
of Aa r is so small that it can be assumed to be negligible and 
Aat is then equal to Aa (fig. 10). 

Waves generated by the addition of heat to adjacent flows of 
different total pressure, considered as case IV (fig. 4), are shown 
in figure 15 at point G- and are illustrated in detail in figure 16(b). 
The strengths of the generated waves, which are equal in magnitude 
but opposite in sign, are obtained from figure 12, depending on the 
value of AQ/Q and Ng • For the case cited, Nq = 10.52, N 2 = 10.39, 
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AQ/Q = 0.03, and Acq is evaluated as -0.0050. In general, the 
value of Aa^ is very small unless Ny and N 2 differ appreci- 
ably. The change in the value of N across the heat-addition zone 
and the waves are obtained successively from figures 5 and 13 for 
the appropriate local values of N, and the corresponding values of 
3 are obtained from figure 14. 

An enlargement of the region near point H in figure 15, in 
which the combination of cases II and IV occurs, is presented in 
figure 16(c). The wave Aa = -0.966 crosses the heat-addition ref- 
erence line at point H, as does the streamline that separates two 
regions of differing total pressure. From the graphical solution of 
case II, the values of the transmitted and the generated waves are 
Aak = -0.995 and Aaj = 0.029, respectively. For case IV, the values 
of the generated waves are Aa^^ = -0.002 and Aa lower = 0.002. 

Therefore, the values of the transmitted and the generated waves of 
case II and case IV are summed algebraically and plotted as -0.997 
and 0.031, respectively, as shown in figure 16(c). 

It is evident that the number of waves and streamlines sepa- 
rating regions in the flow increases as the plot of the flow develops. 
When it becomes desirable to limit the number of waves and stream- 
lines that must be plotted, several devices may be used for combining 
waves and reducing the number of streamlines separating flows of dif- 
fering total pressure. A typical example of the reduction in the 
number of waves by combination of a weak and an adjacent strong wave 
is shown in the region about I in figure 15, wherein weak waves are 
designated by short dashed lines. Figure 17(a) presents an enlarge- 
ment of the portion of figure 15 near point I. The wave Aa = 0.030° 
is combined with the wave Aa = -0.925° by using algebraic addition 
of the wave strengths to give a resulting wave of strength Aa = -0.895°. 
This device is generally employed at the heat-addition reference line 
because every wave that crosses a heat-addition line generates an addi- 
tional wave of minor strength. Figure 17(b), which gives an enlarged 
view of figure 15 near point J, depicts a similar situation where the 
incoming minor wave Aa = 0.034° is combined with the wave Aa = -1.00° 
arising at the wall boundary. The wave that is drawn downstream from 
point J at the wall boundary is equal to -. 09660 } which represents the 
algebraic sum of both wave strengths. This method of combining waves 
avoids the accumulated error that would result from neglecting the 
minor waves and introduces only minor local errors in the flow. 

A method for reducing the number of streamlines separating regions 
of differing total pressures is presented in figures 17(c) and 17(d), 


1387 


NACA TN 2206 


21 


which show an enlarged, portion of figure 15 near point K. In fig- 
ure 17 (c), the streamlines E and F, if drawn past waves U and V, 
would separate regions of different total pressures but of the same 
static pressure. A single streamline H replaces E and F (fig. 17(d)). 
The streamline H is positioned according to the equation 


The flow zone N2 is eliminated. This procedure may somewhat dis- 
tort the rotational properties of the flow; also, when the differences 
in N become large, the accuracy of the resulting plot may be impaired. 

If the number of strong waves that arise at the fixed or free- 
stream boundaries becomes very great, the method illustrated in fig- 
ure 18 can be used at the expense of accuracy in the plot of the flow 
to reduce the plotting work. Along a line immediately upstream of 
a heat-addition zone, a break is made in the plot of the flow. In 
the space provided by the break, a simplified form of the flow is sub- 
stituted for the upstream flow. In figure 18 the letters A, B, C, D 
designate waves immediately upstream of a heat-addition zone. In 
each region of constant total pressure, such as that enclosed by 
streamlines R and S, waves of the same family as A and C are combined 
by adding their strengths. The combined wave is located at the 
weighted mean distance between wave A and C , Aa = Cc, when A and C 
are taken to represent the absolute magnitude of the strengths of the 
corresponding waves. 

The static pressure at any point in the flow field may be obtained 
from the relation 


where P and N are the local values of the total pressure and the 
Mach number squared, respectively. The local value of N appears 
on the graphical plot. Because the total pressure changes with heat 
addition in a manner that depends on the variation of N along each 
streamline, the local value of the total pressure must be obtained by 
tracing the changes in the total pressure along the particular stream- 
line passing through the given point. In order to obtain the change 
in the total pressure across each heat addition of AQ/Q, figure 7 
or the relation 


(N 2 -N 1 )d = (N 3 -N 2 )e 


(27) 


7 
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£2 _ j 1 + 7 w l) 
Pi ~ (1 + 7W2) 



y 

7-1 


may be used, where the subscripts 1 and 2 refer to the regions 
upstream and downstream of a heat-addition reference line on the 
graphical plot of the flow field. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, May 31, 1950. 
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Unheated flow 



Figure 1. - Schematic representation of wave generation at boundary 
between unheated and heated flow; case I. 
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Figure 2. - Schematic representation of wave propagation across heat-addition reference line; case II 
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Schematic representation of wave propagation across adjacent stream 
tubes of differing total pressure; case III. 
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Figure 4. - Schematic representation of wave generation at intersection of heat-addition 
reference line with adjacent stream tubes of differing total pressure; case IV. 
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Figure 6. - Variation of ratio of static-pressure change across heat addition to upstream static pressure with upstream Mach number 

squared. 
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Figure 7. - Variation of ratio of total pressure before heat addition to that after heat addition with upstream Mach number squared. 
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Figure 8. - Variation of strength of wave generated at free-stream boundary of heated flow (from equation (12)) with Mach 

number squared upstream of heat addition; case I. 
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Figure S. - Concluded. Variation of strength of wave generated at free-stream boundary of heated flow (from equation (12)) 

with Mach number squared upstream of heat addition; case I. 
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(a) For heat addition AQ/Q = 0.02 and negative values of incident-wave strength with N]_ = Ng. 


Figure 9. - Variation of strength of wave generated at intersection of incident wave and heat-addition reference line with upstream 

Mach number squared for various incident-wave strengths; case II. 
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(b) For heat addition = 0.02 and positive values of Incident-wave strength with Nq = N2» 

Figure 9. - Continued. Variation of strength of wave generated at Intersection of Incident wave and heat-addltlon reference line 

with upstream Mach number squared for various incident-wave strengths; case II. 
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(c) For heat addition of - 0.03 and negative values of incident-wave strength with N]_ = N 2 « 



Figure 9. - Continued. Variation of strength of wave generated at intersection of incident wave and heat-addition reference line with 

upstream Mach number squared for various incident-wave strengths; case II. 
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(d) For heat addition AQ/Q = 0.03 and positive values of Incident-wave strength with Nq = Ng. 

Figure 9. - Concluded. Variation of strength of wave generated at intersection of incident wave and heat-addition reference line 

with upstream Mach number squared for various incident-wave strengths; case II. 
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Figure 10. - Variation of ratio of transmltted-wave strength to incident-wave strength 
and ratio of ref lected-wave strength to incident-wave strength with ratio A^/Ag, 
where A = YN/VN-1; case III, equation (23b). 
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(a) For heat addition AQ/Q, = 0.02 and values of ^^S.O. 

Figure 12. - Variation of strength of wave generated by heat addition across stream tubes of differing total pressures with Mach 
number squared of region 1 for various values of Mach number squared of region 2; case IV. 
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(b) For heat addition AQ/Q = 0.02 and values of N2 = 3.0. 


Figure 12. - Continued. Variation of strength of wave generated by heat addition across stream tubes of differing total pressures 
with Mach number squared of region 1 for various values of Mach number squared of region 2; case IV. 
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(c) For heat addition of 0.03 and values of NgS'S.O. 

Figure 12. - Continued. Variation of strength of wave generated by heat addition across stream tubes of differing total pressures 
with Mach number squared of region 1 for various values of Mach number squared of region 2; case IV. 
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(d) For heat addition AQ/Q, = 0.03 and values of N£ § 3.0. 


Figure 12. - Concluded. Variation of strength of wave generated by heat addition across stream tubes of differing total pressures with 
Mach number squared of region 1 for various values of Mach number squared of region 2; case IV. 
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Figure 15 - Example shoeing use of method in obtaining graphical solutions for two-dimensional supersonic floes with heat addition. 
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(a) Section about point F, showing wave crossing adjacent 
stream tubes of differing total pressure; case III. 
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(b) Section about point G, showing addition of heat across 
adjacent stream tubes of differing total pressure; case IV. 

Figure 16. - Enlarged sections of graphical example shown in figure 15, 
illustrating in detail several cases. 
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(c) Section about point H showing intersection of wave and streamline separating flows of 
differing total pressure with heat-addition reference line; cases II and IV.- 


Figure 16. - Concluded. Enlarged sections of graphical example shown in figure 15, illustrating in 

detail several cases. 


ui 

co 


NACA TN 2206 


60 


MCA TN 2206 



(a) Section about point I showing method of reducing number of 
waves in flow (combination of minor wave with strong wave). 
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(b) Section about point J showing method of reducing number of waves in flow 
(combination of minor wave with strong wave caused by deflection of wall 
boundary) . 

Figure 17. - Enlarged sections of graphical example shown in figure 15, presenting 
methods for reducing complexity of flow pattern. 
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(c) Section at>out point K before combination of streamlines 
E and F into a single streamline. 



(d) Section about point K after combination of streamlines B and F into 

resultant streamline H. 

Figure 17. - Concluded. Enlarged sections of graphical example shown in figure 15, 
presenting methods for reducing complexity of flow pattern. 
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Figure 18. - Sketch of method of reducing number of waves. 
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